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The analysis of the model quantum clocks proposed by Aharonov et al. Q requires considering 
evanescent components, previously ignored. We also clarify the meaning of the operational time of 
arrival distribution which had been investigated. 
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The concept of time has a very problematic status in the development of quantum theory. This has led several 
researchers to investigate and propose quantum clocks, i. e. quantum devices that in some way will capture and 
reflect a given aspect of time. A recent such set of proposals for quantum clocks, in the specific case of measurement 
C\J ■ of times of arrival, has been put forward by Aharonov ct al. Q|. We shall concentrate here on the first quantum clock 
q \ proposed in that paper, which is coupled to the otherwise free system whose times of arrival one desires to measure 
1} ■ in such a way that the total Hamiltonian (system plus clock) is 

Q: 

<n: h = 2^ p '+ e (- x » p »- d) 

The system variable is x, while the clock corresponds to the time variable y. This is in fact a cyclic variable, which 
entails the conservation of the energy P y . Let us then consider P y restricted to a given value p, in the classical case. 

The motion of the particle (the system) is not fully free: moving from left to right, it runs into a step potential at 
x = 0, either downwards if p > 0, or upwards for negative p. For p > 0, the classical y variable encodes how much time 
elapses since the particle is released and the clock is started, to the instant the particle crosses x = 0. For negative 
barrier height p, the classical particle would be reflected if its energy were not big enough to overcome the step, and 
t-H . the clock variable y would keep running after the particle is reflected. This suggests using predominantly positive p 
values in the analysis, even though quantum mechanically there is reflection even for downward steps. 

In the quantum case, restricting ourselves to an eigenspace of P y with eigenvalue p, either positive or negative, 
the (generalized) eigenstates of the (restricted) Hamiltonian are i) scattering states, with degeneracy two; and ii) 
Qh "evanescent" states, whose eigenvalue is not degenerate. On choosing an adequate orthogonal basis of these scattering 
and evanescent states, one can simply write the time evolution of any given state. However, for states with support 
in the x space restricted to one side (x < 0, in particular), there is a compact alternative expression in terms of an 
integral over a path in complex momentum space which provides us with the whole spacetime dependence of the state 

!!■ 

More explicitly, consider an initial wavefunction in (x,y) space, assumed to be factorized, tp(x,y,0) = ■0i(x) - 02(y), 
. £h ' suc h that ipi(x) has no support on positive x, and compute its Fourier components, that is to say 

' i r + oo p+oc 

^ y '° ) = 2^hJ dk J dpe^ h e^ h g(k)f(p). 

It is then the case that for any time t, the state evolved with the total Hamiltonian H, defined in Eq. (Q), can be 
written as 

/+oo P 
dp j dk f(p)g(k)<pkp(x,y,i), (2) 
-oo Jr(p) 



in 
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where 



, ( ( e ikx/h i k^q_ e -ikx/h\ e ipy/h e -t(p+k 2 /2m)(/S x < 

<t>k P {x,y,t) = — - x { \ k+q ' (3) 

2nn _Zk_ e iqx/h ipy/h -i{p+k /2m)t/h X > 

[_ k+q ' — 



q = \Jk? + 2mp, defined with a branch cut in the k plane, which, for positive p, goes from —i^JImp to i^/2mp, and 
T(p) is a path in the complex k plane from real negative to real positive infinity which goes above the branch cut. 

Thus, in order to write the general solution for Schrodinger's equation with the Hamiltonian H and initial support 
at x < 0, it is imperative to take into account the contribution of the branch cut in Eq. (||), which was omitted in [Q. 
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It should be noted that supposing that g (k) is zero along the branch cuts would in turn bring in the problem that g(k) 
would be forced to be either zero everywhere or non analytic; however, if the initial wavefunction is normalizable and 
has initial support at x < 0, its Fourier transform g(k) must be analytic in the upper half plane. It follows that we 
cannot consistently assume both initial localization and that g(k) is zero along the branch cuts, and the contribution 
of the branch cuts is of necessity present. Quite another issue is whether the contribution of the branch cut can be 
neglected with regard to the physics that we want to describe. If the state were initially peaked at high energies, it 
would present negligible overlap with evanescent states, and their contribution to the posterior evolution of the state 
would remain ignorable. However, we should point out that generically the evanescent component cannot be ignored, 
even in the asymptotic limit t — > oo. As a first example, to be supplemented later graphically, consider the total 
probability that a particle, which starts at t = from the left hand side, x < 0, has to be found at positive x for large 
times, after colliding with a downward step barrier of height p. This (transmission) probability tends to 



dq 



2q 



9(k) 



(4) 



where q and k are related as before, and the evanescent component provides the lower part of the integration interval, 
from to y / 2mp (this expression can be obtained by using a variant of Ricmann-Lebesgue's lemma, in a similar way 
to the computation performed by Allcock in the free case B). Note that there is contribution to (|J) from positive 
and imaginary k, but not from negative values of k. 
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FIG. 1. Collision with a downward step potential: The continuous line portrays the modulus squared of the contribution of 
evanescent waves to the wavefunction; the dashed line is assigned to the modulus squared for scattering waves; the dotted line 
corresponds to the total probability density. All quantities are displayed in atomic units. A time t = 10 a.u. has elapsed. Step 
height p — 2 a.u. . The initial state is a truncated sine (see text), with a — —2.01 a.u. and b — —0.01 a.u. . 



In order to illustrate graphically the relevance of evanescent waves, let us first examine figure |l|, where we depict 
the probability density for positive a; at a given instant, as well as the modulus square of the contribution to the 
wavefunction of evanescent components (corresponding to imaginary k along the branch cut in (^|)), and the modulus 
square of the contribution of scattering (real k) components, in the case of a step potential (fixed p); that is to say, 
we are computing the result of a collision of a wavefunction with a downward step potential. The initial state is 
a truncated sine function, i.e., ipo(x) — (2/(6 — a)) 1 / 2 sin(7r(x — a)/(b — a)) (9(x — a) — 9(x — &)), where 9{x) is the 
step function and > b > a. It is apparent that evanescent waves are numerically important, and even more so the 
interference term between the evanescent and scattering components. 

From classical considerations, Aharonov et al. were led to consider as an operational candidate for the distribution 
of times of arrival of the x particle the following expression: 

/>oo 

Pc (y,t) = / dx\^(x,y,t)\ 2 . (5) 
Jo 
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We would suggest that this quantity be better understood as an operational distribution of dwell times conditional 
on the particle being found at positive x. The corresponding unconditional operational distribution would be given 
by 

/+oo 
dx\ip(x,y,t)\ 2 . (6) 
-oo 

The definition of the unconditional distribution makes apparent its normalization to unity. However, if the evanescent 
components were not included, there would be a probability deficit. Carrying back this argument to the conditional 
distribution (|^), which need not be normalized, one sees that missing out the evanescent components leads to a 
probability deficit. This is apparent in Fig. @. 
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FIG. 2. The conditional distribution (|5|) without evanescent components is represented with a dotted line; with evanescent 
components with the dashed-dotted line. The dashed line corresponds to the unconditional distribution ([]) without evanescent 
components; the continuous line depicts the full unconditional distribution, with the evanescent components taken into account. 
The initial state is a truncated sine wave for the particle, as before, in x space, and a minimum-uncertainty-product Gaussian 
for the clock, with central p being 2 a.u., (y) — 0, and Ay = 1.1 a.u. 



Another effect of considering or not considering the evanescent components is that the peak of arrivals is slightly 
shifted towards higher values of y. This can be understood as due to the fact that the evanescent components are, 
in a way, "slower" than those going over the step, thus remaining longer in x < 0. This also brings down the tail in 
negative y coordinates. 

The additional conceptual distinction between arrival and dwell time distributions is due to the fact that, even in 
a classical picture, particles that arrived at x = but did not cross this point, being reflected, would force the clock 
dial y to keep on moving. Similarly, particles that had crossed x = and were later reflected back would have the 
corresponding dial stop for a while and then resume its ticking. In this way both p c and p u are (conditional and 
unconditional) distributions of dwell-times of the particle in the left half line. 

Another way of understanding the need of reservations with respect to the interpretation of p c as a time-of-arrival 
distribution comes about because of Allcock's analysis of the possibility of ideal distribution of times of arrival (J] . He 
identified the final (t — > oo) probability of finding the particle in the half-space of positive x with the final probability 
of having arrived at x = 0, for a particle initially confined at negative x. He concluded from this identification 
that an ideal time-of-arrival probability distribution could not exist; however, it is by now well known that an ideal 
time-of-arrival distribution does indeed exist for the free particle, namely, Kijowski's distribution [pHj. As has been 
discussed elsewhere [Q]8), the flaw in Allcock's argument was to ignore that negative momenta components (in the 
free particle case and in the case at hand too) would contribute to the arrival probability in a transient regime, thus 
invalidating the stated identification. Similarly, p c (t — ► oo) only contains asymptotic information on dwell times, but 
not detailed information on the transients, so it cannot generically be interpreted as a distribution of times of arrival. 
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FIG. 3. The full conditional distributions at times t = 5 a.u. and t = 7 a.u., for the same initial wavefunction as in Fig. 
(g), are shown with the dotted line. The full unconditional distribution at t = 5 a.u. is depicted with a dashed line, while the 
continuous line represents the full unconditional distribution at t = 7 a.u. . 



The major difference between the conditional and unconditional distributions, when enough time has passed for the 
situation to be considered asymptotic, is the presence in the unconditional distribution of a second peak, associated 
to those components of the wavefunction that are still in x < and will in fact remain there. This second peak is 
present no matter whether we do or do not take into account the evanescent components, and, in fact, does not get 
any contribution from them, as can be seen in Fig. (Q). The evolution of the second peak, moving forward in y space 
as time goes by, is represented in Fig. (||). 

At any rate, it should be noticed that the numerical change of the conditional distribution when the evanescent 
components are taken into account, in the asymptotic regime, is mostly a scale change, thus rendering the difference 
irrelevant for the analysis carried out in with regard to uncertainties in the position of the clock dial. In particular, 
the lower bound suggested in ffl for the product of the average energy of the particle and the uncertainty in the time- 
of-arrival has been demonstrated for Kijowski's distribution without recourse to any coupling with clock variables 
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